DMRG Simulation of the SU(3) AFM Heisenberg Model 
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We analyze the antiferromagnetic SU(3) Heisenberg chain by means of the Density Matrix Renor- 
malization Group (DMRG). The results confirm that the model is critical and the computation of 
its central charge and the scaling dimensions of the first excited states show that the underlying low 
energy conformal field theory is the SU(3)i Wess-Zumino-Novikov-Witten model. 
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I. INTRODUCTION 



In recent years, a renewed interest in models of con- 
densed matter with a symmetry larger than SU(2) has 
arisen. This is because these models represent not only 
challenging theoretical problems but also can be effec- 
tively implemented experimentally. In particular SU(4) 
systems can be realized in laboratories in transition met- 
als oxides [l| where the electron spin is coupled to the or- 
bital degrees of freedom. A possible realization of SU(3) 
antiferromagnetic (AFM) spin chains in systems of ultra- 
cold atoms in optical lattices has been recently proposed 
In this case the spin would be related to the SU(3) 
rotation in an internal space spanned by the three avail- 
able atomic states (colors, in the SU(3) language), with 
the condition that the number of particles of each color 
is conserved. Other examples involve the SU(3) trimer 
state in a spin tetrahedron chain [1, |3| , or the spin tube 
models in a magnetic field Q where the low-energy ef- 
fective Hamiltonian can be identified with a particular 
anisotropic SU(3) spin chain. 

From a theoretical point of view, the SU(3) spin model 
has also been studied from different viewpoints. In recent 
years the interest on ferromagnetic SU(N) spin chains 
has been boosted by their implication in the AdS/CFT 
correspondence [l, |7|. On the other side the family of 
integrable spin chains include some models with SU(3) 
symmetry, as first shown by Sutherland [8], who gener- 
alized the Bethe-Ansatz to multiple component systems 
which include the SU(3) spin chain, showing that it is 
gapless. Also the SU(3) Heisenberg model can be di- 



rectly related to a particular SU(3)-symmetric bilinear 
biquadratic spin-1 chain, the Lai-Sutherland (LS) model, 
which is also known to be critical i^, In terms of 
Conformal Field Theory (CFT) the LS model and the 
SU(3) spin chain should belong to the same universality 
class, that of the SU(3)i Wess-Zumino-Novikov-Witten 
(WZNW) model 

In this paper we present a numerical analysis of the 
SU(3) spin chain by means of the Density Matrix Rcnor- 
malization Group (DMRG). After a short description of 
the model and its mathematical framework (Section [IT]), 
we present our new results (Section IIIip which confirm 
the criticality of the model as well as its correspondence 
to the Lai-Sutherland model. In particular, due to the 
ability of our program to provide the quantum numbers 
for each state, we can show that the excited states of the 
spin chain have the same quantum numbers as the irre- 
ducible representations (IR) of SU(3). We compute the 
scaling dimensions of the first excitations which turn out 
to agree with those of the SU(3)i WZNW model which 
corresponds to the low energy effective field theory de- 
scriptions of our spin chain. The results are further con- 
firmed by the computation of the central charge by means 
of the vacuum entanglement entropy. 



II. THE SU(3) MODEL 

We consider the following Heisenberg model 



H — J Si • Si+i 



(1) 



i=l 
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where the spin variables are expressed in terms of the 
generators of SU(3) in the fundamental representation: 



2 



= iA°, with a = 1,..,8 and Aa the eight Gell-Mann 
matrices. The sign of J selects respectively an antifer- 
romagnetic spin chain (J > 0) or a ferromagnetic (FM) 
one (J < 0). In the following we shall concentrate only 
on the AFM case, which has been partially considered 
also in ref. 

In terms of the following ladder operators, = A"'^±zA^, 
= X^± iX^ and C/± = A^ ± iA^, the Hamiltonian Q 
becomes 



1=1 



4- "'^ I -'^ \8\8 

~^ 2 +^ 2 ' 



(2) 



This makes easier to identify two operators, Sz and Qz, 
given by the sums of the two diagonal Gell-Mann matri- 
ces 



(3) 



that commute with the Hamiltonian and correspond to 
conserved quantities (isospin and hypercharge) . The cor- 
responding quantum numbers label the different eigen- 
states of (dl. 

The Lai-Sutherland model is defined as the bilinear 
biquadratic spin-1 chain 



F = J' 



(4) 



and characterized by an SU(3) symmetry. The model (|4|) 
and the SU(3) spin chain can be mapped one onto the 
other by means of the following identity [isj 



11*^ 

SiSi+i -I- (SiSj+i)^ ~ ^ ^ q + 9 E '^i'^i+i- 



(5) 



a=l 



We have already mentioned in the introduction that the 
LS model is known to be gapless and to belong to the 
same universality class of the SU(3) level-1 Wess-Zumino- 
Novikov-Witten model with central charge c ~ 2. Due to 
the correspondence between the two models, the SU(3)i 
WZNW model has to be the low energy effective criti- 
cal field theory also for the SU(3) spin chain. We shall 
numerically show that the SU(3) Heisenberg chain is crit- 
ical, and from the energy state obtained from the DMRG, 
we shall compute the central charge and the scaling di- 
mensions of (H)) and compare them to the values predicted 
for the SU(3)i WZNW model. 

The states of the spin chain can be organized accord- 
ing to the irreducible representations of the affine (Kac- 
Moody) Lie algebra associated to SU(3). Let us recall 
that a useful way of representing the IR's of the Lie 
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(1) 
(3) 

(3) 
(3) 
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(1) 
(1) 

(1) 

(3) 



(0,0) 

r (±1/2,1/2) 

I (0,-1) 
f (±1/2,-1/2) 
I (0,1) 

(±1/2, ±3/2) 
(0, 0) (3 times) 
(±1,0) 





1/3 
1/3 

2/3 



Table I: Quantum numbers and scaling dimensions for some 
of the primary fields 'I>;,_^ of the SU(3)i WZNW model. 



algebra su(3) is through the Young Tableau (YT) which 
can be labelled by two positive integer numbers (p, g). 
Once p and q are known, one can easily compute the di- 
mension d of the representation and the quantum num- 
bers associated to the isospin Sz and the hypercharge Qz 
according to [H, [l3| : 



d=-{p+l){q + l){p + q + 2) 



and 



Sz 



1,^ 



Qz = lY 



(6) 



(7) 



where / = i(r ± s) and Y = (r — s) — |(p — g), with 
0<r<p,0<s<q. In particular, the cases (1,0) 
and (0, 1) give respectively the fundamental (3) and the 
anti-fundamental (3) IR, while the singlet representation 
(1) corresponds to (0^ 0). 

It has been proved [18|] that, in analogy with the SU(2) 
case, the ground state (GS) of the AFM SU(3) Hamil- 
tonian is a singlet and, since it is made of particles m, 
d and s in equal number, it can be obtained in finite 
chains having only a number of sites which is a multiple 
of three, L = 3M. As for the excited states, we expect 
them to be in correspondence with the tower of conformal 
states of the corresponding SU(3) WZNW model. The 
primary states of this theory are a finite number and are 
given [l^l by fields whose holomorphic (antiholo- 

morphic) part transforms according to a representation 
A = (p, q) (A = (p', q)) with the values of p, q (and simi- 
larly of p', q') satisfying the condition: p + q < k, k being 
the level. The conformal dimension of the primary field 



is then ^ 



with 



i{k + N) 



{p'' + q^+pq + 3p + 3q), (8) 



and a similar expression for x 



{p',q') 



For future reference, 



the values of x^ \ for some primary fields in the case of 
k ~ \ arc reported in Table HI 

To end up this section, we notice that in a finite chain 
of length L not all quantum numbers, i.e. states, may be 
realized. For examples, working with periodic boundary 



conditions and with an even number of sites, the singlet 
(1) X (1) (ground) state, with x — 0, appears only for 
chains with L = 6M (with M a positive and integer 
number), while the (3) x (1) (or the (1) x (3)) states are 
present only if L = 6M + 4 (or L = 6M + 2), both with 
X = 1/3. 



III. NUMERICAL ANALYSIS 

The SU(3) version of the DMRG we have used imple- 
ments the following Hamiltonian 



H 



T ^ 



(9) 



+ k2 U+U^+^ + h.c.) + - {zo XfX^+i + zi 



where kj and 



are input parameters. The model ^ 
reproduces the AFM (FM) case when all the fcj's and 
the 2j's are equal to 1 (-1). By tuning the input param- 
eters kj and Zj , we can study all the possible anisotropic 
version of the SU(3) Heisenberg model. A very impor- 
tant feature of this DMRG is that it implements both 
the quantum numbers Sz and Qz given in ([3]). This im- 
plementation considerably reduces the computation time 
and, on the other hand, once Sz and Qz are fixed from 
input, each run of the DMRG yields exclusively the ener- 
gies of the states within those quantum-number sectors. 
This is very useful when one needs to classify the exci- 
tations according to the values of the isospin and of the 
hypercharge. 

By setting ki — k2 = Zi = 0, we restrict to the SU(2) 
sector of SU(3). This has been used as a check to the 
program; the DMRG in this case reproduces perfectly all 
the energy states of the SU(2) Heisenberg model. 

We study now the isotropic AFM chain with periodic 
boundary conditions by means of an infinite size DMRG 
with up to TO = 2200 states in order to reduce the uncer- 
tainty on the energies to the order of magnitude of the 
truncation error. The data for the ground state and the 
first excited states are plotted in Fig. [1] 

Let us first concentrate on the ground state, which, in 
agreement with theoretical predictions, it is found only 
when L = 6M. The plot of Eqq as a function of 
shows a good linear behavior; this justifies the fitting of 
our data by the CFT equations for the GS: 



L 



ncv 



(10) 



where Coo and the product cv are kept as fitting param- 
eters. We obtain: Coo = -0.518288 and cv = 2.04419. 
In order to derive the value of v we need an independent 
derivation of c. The central charge for a SU(N) level-k 
WZNW model is given by [11] 



k {N^ - 1) 
k + N 



(11) 
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Figure 1: Plot of the ground state Eoo and of the excited 
states for chains of different lengths (from L = 12 up to L = 
52). The solid lines have a slope respectively of i and | and 
have been drawn as a guide for the eye. For sake of clarity, 
not all the degeneracies have been reported. 



If the effective field theory describing our spin chain is 
the conformal SU(3)i WZNW model, the central charge 
must be c = 2. 

However, it is possible to have a direct numerical 
derivation of c from the asymptotic behavior of the von 
Neumann entropy Sn = — Tr„(p„ log2 Pn) of the reduced 
density matrix p„ = Tri>„p of a subchain with n spins 
of a critical system of length L, as a function of n and L, 
where p is the density matrix associated to the ground 
state of the chain. Indeed, one has [19, 20]: 



■l0g2 



L . fTT \ 

— sm —n 

TT \L J 



A. 



(12) 



As usual c is the central charge while A is a non-universal 
constant. The DMRG computes the density matrix for 
a block of length n in a chain of length L, so that Sn 
becomes quite simple to calculate. Fig. [5] shows the be- 
havior of the von Neumann entropy as a function of the 
block length n in (a) and as a function of the quantity 
y — log2 (^) /3 (obtained from by setting n — L/2) 
in (b), for values of the DMRG states equal to 1000 and 
2000. The data confirm the linear behavior expected 
from Eq. (fT2|) . The linear regression Sn = cy ~\- A yields 
the value for the constant A = (1.774±0.002) and for the 
central charge c = (1.995 ± 0.001). Thus the theoretical 
prediction of Eq. (jlip is confirmed with very high accu- 
racy. From Fig. [H^b) it is also evident that the values 
obtained when keeping only 1000 states in the DMRG 
run are much less precise. This is the reason why we 
have then performed all calculations while keeping 2000 
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0.5 1 
log2(L/7i)/3 



L 




X2 


6M+2 
6M+4 
(A) 


0.3414 ± 0.0001 
0.3406 ± 0.0002 
0.3410 ± 0.0002 


0.6291 ±0.0003 
0.6238 ± 0.0003 
0.6265 ± 0.0003 


6M 




0.6503 ± 0.0003 



Table II: Results of the numerical analysis. The velocity and 
the central charge are respectively: v = (1.0247± 0.0005) and 
c = (1.995 ± 0.001). The scaling dimensions of the first (a;i) 
and the second (3:2) excited states for each L are calculated 
as described in the text. The mean value (a) between L = 
6M + 2 and L = 6M+4 for x\ and X2 is provided (see also Fig. 
[T]). For L = 6A/ only the first excitation above the ground 
state has been considered. 



Figure 2: Analysis of the von Neumann entropy. The figure 
on the left (a) shows Sn for m = 1000 DMRG states as a 
function of the block length n for different chains of length L; 
the figure on the right (b) is the plot of the half chain entropy 
(n = L/2) for m = 1000 and m = 2000 DMRG states. The 
linear fit on these data provides the values for the central 
charge c and the constant A. 

states. Finally, the value of c can be substituted into the 
product cv derived from the GS to recover the velocity of 
the excited modes: v = (1, 0247 ± 0.0005), which is close 
to the expected value [8| 7r/3. 

Before proceeding with the analysis of the excited 
states, let us check the asymptotic value of the energy 
density Eoo- The theoretical prediction for the ground 
state of the S = 1 bilinear biquadratic Heisenberg Hamil- 
tonian (see Ref. [2l[) is 

i?G5 = -ln3-^ + l, (13) 

which already takes into account the factor -1 of the 
l.h.s. of equation Starting from the correspondence 
between our SU(3) chain and the biquadratic one (O, 
we can compare the value of we obtained with the 
one predicted by equation Eqs = -0.703212. The 

match is exact to the third decimal (—0.703243), if one 
also recalls that the Hamiltonian has a factor 1 /4 (due to 
the definition of the spin variables in terms of the SU(3) 
generators) so that Coo needs to be multiplied by a factor 
two, and summed to the factor 1/3 of equation ([5]). This 
is a further numerical proof of the equivalence between 
the the Lai-Sutherland and the SU(3) spin model. 

Let us study now the excited states. From Fig. [U one 
immediately sees that the slope of excited states depends 
on the length of the chain. In particular, for L = 6M the 
first excitation scales with a slope which is unmistakably 
different from the slope of the L = 6M + 2 or L = 6M + 4 
data. For small values of L the data corresponding to the 
same Sz but with opposite Qz are split by a finite size 
correction, while for increasing values of L they tend to 
overlap and scale to the same asymptotic value. 

For the excited states CFT predicts that: 

E, - Eoo = (14) 



where Xj is the scaling dimension of the j— th excitation 
for a given chain of length L; Eqq is given by Eq. (fTO|) 
where c and v have been derived before and are reported 
in the caption of Table [ill The numerical coefficients 
for the scaling dimensions that one can obtain from the 
DMRG data of Fig. [T]are hsted in Table HH 

As expected, the values of the allowed conformal di- 
mensions are very close to the values of 1 /3 and 2 /3 pre- 
dicted by a SU(3)i WZNW model. 



IV. CONCLUSIONS 

We have provided strong numerical evidence of the 
criticality of the AFM SU(3) spin chain. Also, we have 
confirmed that the conformal field theory describing the 
chain is effectively the SU(3)i WZNW model, by compu- 
tating the central charge and scaling dimensions of the 
lowest excited states of the model, which turn out to be 
organized according to the IR of SU(3)i Kac-Moody al- 
gebra. 

There are many interesting generalizations of the above 
models which deserve further study. In particular, a sim- 
ilar ferromagnetic spin chain is connected with the non- 
linear C'P^ sigma mode at 6 = ir and might be useful 
to clarify some controversial problems of the model. An- 
other interesting problem is to consider larger SU{N) 
symmetry groups. In two-dimensional chains, the vac- 
uum state is of Neel type for N < 4 and of Spin-Peierls 
type for iV > 5 [13]. The analysis by means of DMRG 
technique might shed some light on the transition mech- 
anism. 
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